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Abstract: Genetic variability in coypus, Myocastor coypus, and its relation with hunting pressure. 

Genetic variability of coypu in sites that support different hunting pressures was analysed through PCR 
amplification of RAPD markers. DNA was extracted from tissue samples obtalned in tive freshwater sites located ¡n 
Luján zone, north-east of Buenos Aires province, Argentina, between October 1999 and September 2000. Fifty-one 
RAPD markers were obtalned. Forty-nine of these were polymorphic between individuáis. The analysis of 
molecular variance (AMOVA) showed a high genetic variability between sites, which would Indícate that individuáis 
from these sites are genetically differentiated. Genetic diversity within each site was negatively correlated with 
hunting pressure, while a positive correlation between genetic diversity and population size was observed. Genetic 
distance between ¡mmigrants, residents and immigrants-residents were significantly different, showing that 
immigrant individuáis are more related between them than with the rest. It is proposed that in the Pampas Región, 
where the landscape ¡s fragmented into pareéis having agricultural-cattle use, a source-sink dynamic between 
pareéis protected and not protected from hunting pressure is establlshed for the coypu where protected pareéis 
would act as a source of individuáis for the not protected ones. 
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Muchas poblaciones de mamíferos se encuen¬ 
tran amenazadas de extinción o han desapareci-do 
a causa de la intensa explotación realizada por el 
hombre (Mace & Reynolds, 2001). La activi-dad de 
caza representa la segunda amenaza a las especies 
de mamíferos, después de la degradación del 
hábitat (Mace & Reynolds, 2001). Sin embar-go, la 
explotación de fauna también representa una 
importante fuente de alimento y de ingreso 
económico para pobladores rurales, particular¬ 
mente en el Neotrópico (Robinson & Redford, 1991). 
El coipo o nutria Myocastor coypus es un roedor 
grande que se distribuye en ambientes acuá-ticos 
del sur de Sudamérica. En Argentina, el coipo es un 
importante recurso económico, ocu-pando el primer 
lugar en las exportaciones lega-les de fauna 
silvestre (Marchetti & Morello, 1991; Bárbaro, 1994) 
. Sólo una pequeña cantidad de pieles proviene de 
criaderos, mientras que el 98% de las pieles 
comercializadas se extraen de pobla-ciones 
silvestres (Zacagnini & Venturino, 1993). Este hecho 
ha llevado a la sobreexplotación de poblaciones de 
coipo en algunas partes de su ran-go nativo 
(Gosling & Baker, 1991), por lo que cons¬ 


tituye una de las especies que requiere de 
progra-mas de manejo para lograr una 
explotación sus-tentable en Argentina (Holdo & 
García Fernández, 1995). 

La perdurabilidad de las poblaciones bajo ex¬ 
plotación depende del efecto de las estrategias 
de manejo sobre la demografía de la especie 
(Taylor & Dunstone, 1996; Boyce et al., 1999). El 
efecto de la caza sobre la variabilidad genética 
de las poblaciones, mediado por la disminución 
del tamaño poblacional, es uno de los principales 
aspectos demográficos que puede afectar la via¬ 
bilidad de poblaciones silvestres. La disminución 
de la diversidad genética puede deberse a diver¬ 
sos factores como endogamia, deriva genética y 
cuellos de botella, cuyos efectos se incrementan 
cuando el tamaño poblacional disminuye 
(Frankham, 1996). Flasta el momento, no se co¬ 
nocen estudios donde se haya estudiado la 
estruc-tura genética de poblaciones o grupos 
sociales de coipo. 

Los objetivos de este estudio fueron (1) poner 
a punto la técnica RAPD (Random Amplified 
Polymorphic DNA) (Williams, 1990) para identi- 
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ficar marcadores para el coipo y (2) analizar la 
estructura genética en distintos sitios utilizados 
por el coipo que reciben presiones de caza dife¬ 
rentes. La técnica RAPD se basa en la amplifica¬ 
ción al azar de fragmentos de ADN por PCR 
(iPolymerase Chain Reaction). Esta técnica per¬ 
mite el procesamiento de una gran cantidad de 
muestras con un costo bajo de tiempo y dinero y 
requiere solo una pequeña cantidad de ADN. 
Debido a que el nivel de polimorfismo detectado 
por los marcadores RAPD es frecuentemente 
ele-vado, estos son útiles en la descripción de la 
es-tructura genética (Bowditch, 1993). 

Se trabajó en cinco sitios, dos con presión de 
caza baja, dos con presión de caza alta y un quin-to 
con una presión intermedia. La hipótesis que se 
puso a prueba fue que los sitios bajo presión de 
caza debían mostrar una menor variabilidad 
genética por el efecto de la caza sobre los núme-ros 
poblacionales, a menos que una elevada tasa de 
inmigración aumentara la diversidad genética. 

MATERIALES Y METODOS 

Sitios de estudio 

El trabajo se llevó a cabo en cinco cuerpos de 
agua ubicados en la zona de Luján, al nordeste 
de la provincia de Buenos Aires, Argentina, en¬ 
tre octubre 1999 y septiembre 2000. Esta es una 
zona típica de la región de la Pampa Ondulada 
en la que la vegetación natural, caracterizada por 
la estepa graminosa, ha sufrido importantes mo¬ 
dificaciones y fue reemplazada por especies exó¬ 
ticas sembradas como cultivos, pasturas o 
forestaciones o que crecen de manera espontá¬ 
nea (Cabrera & Zardini, 1993). El clima es tem¬ 
plado y húmedo, con una temperatura media de 
9.1 °C en invierno y 23.8 °C en verano y una pre¬ 
cipitación media anual de 944 mm (Goldberg et 
al., 1995). 

El primer sitio de estudio, Ruiz (34° 12’ S, 59° 
16’ O), se encuentra en un segmento de arroyo de 
1,5 km de largo. El segundo sitio, Luján (34° 35’ S, 
59° 04’ O), esta formado por dos pequeños 
estanques (0,6 y 0,08 has) separados por una por¬ 
ción de arroyo de 100 m de largo. Los otros tres 
sitios, Jáuregui, Ratti-1 y Ratti -2 (34°35’ S, 59° IT 
O) son tres lagunas (4,3; 2,3 y 4 has, respecti¬ 
vamente) resultantes de la modificación de cur-sos 
de agua naturales. Jáuregui se encuentra ubicada a 
4 km de los sitios Ratti-1 y Ratti-2, que están a su 
vez separados por 200 m. Las tres lagunas están 
conectadas con el río Luján por arroyos o canales 
de menos de 100 m de longitud en el caso de Ratti- 
1 y Ratti-2, y de 500 m de lon-gitud en el caso de 
Jáuregui. 


Existió una diferente presión de caza entre 
sitios, siendo elevada en Ratti-1 y Ratti-2, inter¬ 
media en Luján y nula-baja en Ruiz y Jáuregui. 
Estos resultados fueron obtenidos durante las 
campañas de captura viva y en una o dos visitas 
a cada sitio entre campañas mediante la 
medición en cada sitio de (1) el número de 
madrigueras cavadas utilizando una pala, (2) la 
presencia de cazadores, y (3) signos de actividad 
de cazadores mediante el conteo de cartuchos 
vacíos (Guichón & Cassini, 2005). 

Recolección de muestras 

Se realizaron campañas de captura viva con 
frecuencia estacional en la primavera (Oct -Nov) de 
1999 y en el verano (Ene-Mar), el otoño (Abr-Jun) y 
el invierno (Jul-Sep) de 2000 en todos los sitios de 
estudio (métodos descriptos en Guichón et al., 
2003a; Guichón et al., 2003b). Las sesiones de 
captura tuvieron una duración de entre 4 y 11 días y 
fueron diseñadas para capturar todos los coipos 
presentes en cada sitio, finalizando la cam-paña 
cuando la proporción de recapturas exce-día el 80% 
del total de capturas por dos días con-secutivos . 
Las trampas, ubicadas en la abertura de las cuevas 
o en senderos o sitios de forrajeo cercanos a estas 
fueron cebadas con manzana antes del atardecer, y 
chequeadas antes del ama-necer. Las 
observaciones entre sesiones de cap-tura indicaron 
que la gran mayoría de los coipos residentes en 
cada sitio eran marcados bajo este régimen 
(Guichón et al., 2003b). Los coipos fue-ron 
inmovilizados con una inyección intramus-cular que 
contenía una combinación de clorhi-drato de 
ketamina (4 mg/kg) y clorhidrato de xylazina (0.5 
mg/kg) (Bó et al., 1994) y se les tomo una pequeña 
muestra de tejido de la oreja mien-tras se procedía 
a su marcado mediante carava-nas numeradas de 
distinto color y forma. Los in-dividuos eran 
clasificados como inmigrantes o residentes sobre la 
base de su historia de captu-ra en el sitio y el peso 
en la primer captura. En la primer campaña de 
captura todos ios individuos fueron considerados 
residentes. En las siguien-tes campañas, ios 
individuos capturados por pri-mera vez con un peso 
corporal < 1,5 kg fueron considerados como nacidos 
en el lugar entre cam-pañas de captura, 
basándonos en Brown (1975: los juveniles pesaron 
< 1,25 kg) y nuestros da-tos de recaptura (los 
juveniles que pesaron 0,2-1,5 kg en la primer 
captura pesaron 1,1-2,5 kg cuando eran 
recapturados, mientras que no se capturo ninguna 
cría nueva de 1,5-2 kg). Los demás individuos 
capturados por primera vez fueron considerados 
inmigrantes (Guichón et al., 2003b). 
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Métodos moleculares 

Se seleccionaron muestras de tejido de 24 in¬ 
dividuos de los cinco sitios: ocho de Ruiz (n = 8), 
dos de Luján (n = 2), siete de Jáuregui (n = 7), tres 
de Ratti-1 (n = 3) y cuatro de Ratti-2 (n = 4). La 
extracción del ADN se inició disgregando el tejido y 
colocándolo en un tubo eppendorf con 500 «I de 
buffer de extracción (10 mM Tris-CI, 0.1 M EDTA, 20 
ocg/ml de RNAsa, 0.5% SDS, 120 ocg/ml de 
proteinaza k, pH 8.0). Posteriormente se incu-bó la 
suspensión en un baño a 50°C por 16 hs. 
(Sambrook et al., 1989). Al día siguiente se com¬ 
pletó la extracción utilizando el método de fenol- 
cloroformo (Sambrook et al. , 1989). Posteriormen-te 
el ADN fue cuantificado en espectrofotómetro 
UV/VIS (Hitachi, modelo U-2001) a 260/280 nm. 
Durante la puesta a punto de la técnica se utiliza-ron 
10 cebadores (serie A01 a al 0) (Biodynamics 
S.R.L., Argentina) y se probaron distintas concen¬ 
traciones de todos los reactivos necesarios para la 
PCR para lograr la puesta a punto de la técnica. La 
mezcla de reacción (volumen final 25 ocl) con-sistió 
en: 35 ng de ADN, 100 ocm de cada dNTP, 200 ocm 
de cebador, buffer de reacción IX, 3 mM de MgCI 2 y 
1,5 unidades de Taq ADN polimerasa (Promega). 
Cada reacción fue cubierta con 30 °d de aceite 
mineral para prevenir la evaporación. Para la 
amplificación, las muestras fueron some-tidas a 1 
ciclo de: 2,5 minutos a 94° C, 1 minuto a 35° C y 2 
minutos a 72° C, seguido de 44 ciclos de: 1 minuto a 
94° C, 1 minuto a 35° C y 2 minutos a 72° C en un 
termociclador Perkin Elmer Cetus (Emeryville, Calif). 
Los fragmentos generados por la amplificación 
fueron separados en geles de agarosa al 1%, 
teñidos con Bromuro de Etidio y visualizados y 
fotografiados bajo luz UV. Para po-der comparar los 
patrones de bandas obtenidos en los distintos geles 
y poner a prueba la replicabilidad de los resultados, 
seis de las vein-ticuatro muestras analizadas (el 
25%) fueron amplificadas y corridas 
electroforéticamente por duplicado, obte -niéndose 
siempre el mismo pa-trón de bandas para cada una 
de las seis mues-tras, lo que demuestra una 
replicabilidad de los resultados del 100 %. 


Análisis de datos 

Los datos de 51 marcadores, fueron incorpo¬ 
rados en una matriz de ceros (ausencia de la 
ban-da) y unos (presencia de la banda) en el gel 
de agarosa. Las bandas de igual tamaño fueron 
in-terpretadas como homologas (Gugerli et al., 
1999). Para estimar la variabilidad genética en¬ 
tre los sitios de estudio se realizó un análisis de 
varianza molecular (AMOVA) (Excoffier et al., 
1992), utilizando el programa Arlequín 2.000 
(Schneider et al., 2000). A partir de la matriz de 


ceros y unos se calculó una matriz de distancia 
genética de Nei (Nei, 1972), utilizando el pro¬ 
grama Rapdistance 1.0 (Armstrong, 1996). Esta 
matriz fue la base para los análisis subsecuen¬ 
tes. 

Para comparar la distancia genética entre in¬ 
dividuos inmigrantes e individuos residentes, se 
obtuvieron valores de distancia genética entre: 

(1) pares de inmigrantes, (2) pares de residentes 
y (3) pares de un inmigrante y un residente. Es¬ 
tas comparaciones incluyeron 5, 17 y 27 valores 
de distancia genética, respectivamente. Para 
este análisis solo se utilizaron los datos de 
individuos de los sitios Jáuregui, Ratti-2 y Ruiz, 
debido a que las dos muestras de Luján eran 
sólo de indi-viduos residentes y en las tres 
muestras de Ratti-1 sólo había un inmigrante. La 
distancia genética dentro de cada sitio fue 
estimada mediante el índice de diversidad de 
Shannon-Wiener (Lewontin, 1973): Hs = -Zp¡ ,log 2 
p¡, donde p, es la frecuencia de cada marcador 
RAPD en cada sitio de muestreo (Dawson et al., 
1995). Luego se realizó un análisis de correlación 
de Spearman entre los valores del índice de 
diversidad y pre-sión de caza y entre los valores 
del índice de di-versidad y el tamaño poblacional 
en cada sitio. Para esto se utilizaron los datos de 
tamaño poblacional y presión de caza estimados 
por Guichón & Cassini (2005). 

RESULTADOS 

Se obtuvo amplificación con seis de los diez 
cebadores utilizados (A01, A03, A06, A08, A09 y 
Al 0), dando como resultado 51 marcadores 
RAPD de un tamaño entre 200 y 1500 pb. Cua¬ 
renta y nueve de estos marcadores resultaron 
polimórficos entre todos los individuos (Fig. 1). 

El análisis de varianza molecular de estos 
marcadores reveló una elevada variabilidad 
genética entre sitios (AMOVA, 0 S t = 0,39, n = 
24; P = 0,001), lo que indicaría que los indivi¬ 
duos de estos sitios se encuentran 
genéticamente diferenciados. 

La diversidad genética estuvo negativamen-te 
correlacionada con la presión de caza 
(Spearman, r s = -0,90, n = 5; P = 0,023), mien¬ 
tras que se observó una correlación positiva en¬ 
tre diversidad genética y tamaño poblacional 
(Spearman, r s = 0,97, n = 5; P = 0,012). Por otro 
lado, no se obtuvo una correlación significativa 
entre los valores del índice de diversidad y el ta¬ 
maño de muestra utilizado (Spearman, r s = 0,60, 
n = 5; P = 0,28), lo que indica que los resultados 
anteriores no son la resultante de un artefacto en 
el número de individuos analizados de cada sitio. 
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Fig. 1. Gel de agarosa 1% mostrando las bandas 
obtenidas mediante la amplificación por PCR del 
ADN de los 24 individuos analizados utilizando el 
cebador A10. 


La distancia genética promedio entre 
inmigrantes fue de 0,10 ± 0,04 (n = 5); entre re¬ 
sidentes de 0,22 ± 0,07 (n = 17) y entre 
inmigrantes y residentes de 0,26 ± 0,10 (n = 27). 
Estas distancias fueron significativamente dife¬ 
rentes entre sí (prueba de Kruskal-Wallis, H = 
10,43; P = 0,005). 


DISCUSION 

En el presente estudio se utilizó la técnica de 
marcadores RAPD con el objetivo de analizar la 
estructura genética del coipo en distintos sitios. 
El análisis de varianza molecular indico que la 
variabilidad entre los sitios estudiados es grande 
(AMOVA, 0 ST = 0,39, n = 24; P = 0,001). Esta 
proporción de variabilidad es similar a la encon¬ 
trada en otro roedor, Peromyscus maniculatus, 
para el cual aproximadamente el 36% de la di¬ 
versidad fue encontrada entre poblaciones, utili¬ 
zando la misma técnica (Vucetich etal., 2001). 

Lynch & Milligan (1994) han propuesto que el 
sesgo en la estimación de parámetros 


poblacionales mediante el uso de marcadores 
RAPD puede ser despreciable si se restringe el 
análisis a marcadores con una frecuencia 
observa-da menor a 1 - 3/N, donde N es el 
número de individuos analizados. Sin embargo, 
en dos estu-dios recientes (Bartish et al., 1999; 
Vucetich et al., 2001) no se encontraron 
diferencias en los resultados obtenidos utilizando 
todo el conjunto de marcadores o solo aquellos 
que cumplían con el criterio propuesto por Lynch 
y Milligan (1994). Bartish etal. (1999) concluyen 
que cuando se tra-baja con pequeños tamaños 
muéstrales, como ocurre en nuestro caso, la 
utilización del criterio de Lynch & Milligan (1994) 
podría llevar a la subestimación de la diversidad 
genética. Tenien-do en cuenta esto último 
podemos decir que la técnica aplicada en este 
estudio es un método confiable para analizar la 
estructura genética, si bien está limitada por el 
pequeño tamaño de la muestra analizada. 

Los modelos clásicos de poblaciones bajo ex¬ 
plotación, consideran que los animales extraídos 
son individuos resultantes de la tasa de repro¬ 
ducción dentro del área donde se realiza la caza y 
no tienen en cuenta que estos podrían ser indi¬ 
viduos provenientes de otras poblaciones (Caughley 
& Sinclair, 1994). Modelos más recien-tes 
incorporan la estructura espacial de las po¬ 
blaciones y enfatizan el rol de la dispersión en-tre 
sitios con diferente presión de caza, como por 
ejemplo mediante la aplicación de modelos fuen-te- 
sumidero a poblaciones bajo explotación (Tuck & 
Possingham, 1994; Lundberg & Jonzén, 1999; 
Novara et al., 1999; Milner -Gulland, 2001) . En 
nuestros sitios de estudio, a pesar de que la caza 
no afectó la tasa de reproducción de los coipos, la 
mortalidad en los sitios con caza fue tan elevada 
que la mayoría del reclutamiento se debió a la 
inmigración (Guichón & Cassini, 2005). Es pro-bable 
que bajo este contexto se haya establecido una 
dinámica de tipo fuente-sumidero entre si-tios 
protegidos y no protegidos de la caza. Con¬ 
trariamente a lo esperado, en este estudio la in¬ 
migración no trajo aparejado un aumento en la 
diversidad genética. Esto pudo deberse a un efec-to 
no deseado del tamaño muestral o que los 
inmigrantes provienen de la misma población y 
quizás de la misma sub-población. En los coipos, la 
dispersión está asociada a los cursos de agua por lo 
que es posible que los inmigrantes lleguen desde la 
población más cercana conectada por ríos y 
arroyos. Adicionalmente, los individuos dispersantes 
podrían ser parientes, como pare-ce ocurrir en otras 
especies de roedores cavio-morfos (Herrera & 
Macdonald, 1987; Branch et al., 1993). 
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